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ABSTRACT 

Hybrid control systems are becoming increasingly important today due 
primarily to the development of the high speed digital computer. This 
work involves the hybrid simulation of a relatively simple second order 
analog system using a short word length digital computer for the computa- 
tion of the control law for the analog system. Particular attention is 
paid to the effects of the analog-digital interface on the analog system 
response. A method is developed whereby the digital computer can be 
programmed to minimize these interface effects. 

The simulation was Girried out in the Digital Control Laboratory at 
the U. S. Naval Postgraduate School using the CDC 160 digital computer 
and the Pace TR-20 analog computer with associated equipment. Although 
the specific solution of the interface problem is directed towards the 
above equipment, the general method of attack on the interface problem 


is applicable.to many similar type hybrid systems. 
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CHAPTER I 
INTRODUCTION 

The combined use of analog and digital computers in the simulation 
of complex systems wal increased manyfold in recent years, due primarily 
to the increase in speed of computations of the digital computers. These 
hybrid simulation schemes seein the parallel computation features of the 
analog computer with the accuracy and versatility of the digital computer. 

This study was undertaken with the purpose of developing a series of 
programs and subroutines which could be used for the simulation of a large 
number of hybrid control systems using a digital computer with a short word 
length and limited arithmetic computational ability. As the study develop- 
ed, it became apparent that one of the more important considerations was 
the analog-digital interface and its possible effect on system response. 

As a result, the major effort was directed towards the detailed study of 
a simple second order hybrid system simulation in order to investigate 
these effects in terms of a specific system. 

Chapter two gives a brief description of the equipment used in the 
simulation study; Chapter three describes the general considerations which 
must be considered in any hybrid control scheme, and Chapter four deals with 
the actual simulation study itself. 

Towards the end of developing general programs, a subroutine (Appen- 
dix II) was developed for general use which provides for the introduction 
of computational constants into the digital computer in a format compat- 
ible with the analog-digital conversion equipment. Also, the digital com- 
puter program developed for the system simulation is applicable to any 
second-order analog system where the digital computer is used to compute ‘a 


control value consisting of a linear combination of the states of the 
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analog system. Also, the general method of attack on the limitations 
imposed by the A/D interface is applicable to hybrid control systems in 
general although the specific procedures may vary with each specific 


system. 





CHAPTER 2 
DESCRIPTION OF EQUIPMENT 
2.0 GENERAL. 

The experimental phases of this thesis were carried out in the 
Digital Control Laboratory, U. S. Naval Postgraduate School. System 
simulation was accomplished using the following equipment: 

A. The Control Data Corporation 160 digital computer. 
B. The Control Data Corporation 168 arithmetic unit. 
C. The Pace TR-20 analog computer. 
D. A 12-bit, multi-channel analog-to-digital and digital-to- 
analog conversion unit (ADC-DAC). 
Brief descriptions of the above equipment are given in the following para- 


graphs. 


2.1 CONTROL DATA CORPORATION 160 DIGITAL COMPUTER. 

The CDC 160 digital computer is a electronic computer controlled by 
an internally stored program in sequential locations. Memory capacity is 
4096 12-bit binary words of magnetic core storage, with a storage cycle 
time of 6.4 microseconds. Instructions are executed in one to four stor- 
age or memory cycles with the time required for execution varying from 
6.4 to 25.6 mi¢roseconds. The average instruction execution time is 15 
microseconds. . 

A CDC 160 computer word is made up of 12 binary digits. The bits 


within a computer word are numbered from 0 to 11, right to left, i.e., 


11 10 09 08 07 06 05 04 03 02 01 O00 


All arithmetic is binary, one's complement notation. Although the computer 
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operates in the binary system, the octal representation of a binary num- 
ber is more convenient. The CDC 160 word can then be expressed as four 


octal digits, i.e., 






Sonmee 


All positive numbers must have a "0" in bit 11; all negative numbers 
must have a "1"' in bit 11. In octal notation, numbers from (0000), thru 
(3777), are positive; numbers from (4000), thru (7777), are negative. 

All programs were compiled using the One Sixty Assembly Program (OSAP). 
The OSAP output consists of a side-by-side listing of the machine language 
instruction code and a alpha-numeric program listing. The CDC machine 
language instruction word is divided into a six bit function code (F code) 
followed by a six bit executroni address (E field). Most instructions follow 
this 12-bit pattern and thus require only one word of storage; however, 
certain expanded instructions require 24-bits thus using two words of 
storage. The first word of such a two word instruction consists of the 
six-bit function code followed by a six-bit execution address which is 
always zero. The succeeding 12-bit word contains the address or operand 
(G code) depending on the instruction. 

The CDC 160 computer communicates with external or peripheral equip- 
ment by means of external function codes. A detailed description of these 
external function codes is given in Appendix II. 

The CDC 160 computer is very limited in its arithmetic computation 
capability. Twelve bit addition and subtraction is accomplished in two or 
three memory cycles, depending upon the address mode used. Multiplication 


and division can only be accomplished by successive 12-bit addition or 
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subtraction, respectively, resulting in excessive programming and execu- 


tion time requirements. 


2.2 CONTROL DATA CORPORATION 168 ARITHMETIC UNIT. 
The peripheral CDC 168 arithmetic unit provides the CDC 160 computer 
with the capability of double precision addition and subtraction plus 
single or double precision multiplication and division of data. The CDC 
168 is addressed by the CDC 160 with external function codes (See Appendix 
II) which select the unit and specify which operation is to be performed. 
The computer then transmits the proper number of words in the correct 
format to make up the arguments for the operation selected. The CDC 168 
proceeds to compute the result and awaits the computer request for results. 
In the programs which follow, the CDC 168 is used in the single pre- 
cision multiply and divide modes. Single precision (short) multiplication 
uses an 11-bit multiplier and an 11-bit multiplicand, providing a 22-bit 
product. Single precision (short) division requires a 22-bit dividend and 
a ll-bit divisor yielding a 1l-bit quotient. This operation requires from 
145 to 180 microseconds while the short multiply operation requires 120 
microseconds, The 22-bit arguments are formed from two 12-bit CDC 160 
words using the least significant 1l-bits of the two words while the most 
significant bit in each word is considered as a sign bit. Appendix I gives 
several examples of the short multiply and divide operations. Double pre- 
cision multiplication (22 bits x 22 bits) and division (44 bits + 22 bits) 


are available at the expense of doubling the time per operation. 


2.3 ANALOG-TO-DIGITAL AND DIGITAL-TO-ANALOG CONVERSION 
EQUIPMENT, 
A 12-bit, multi-channel analog-to-digital and digital-to-analog con- 


version package was used to convert analog voltages to digital numbers and 
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to convert digital numbers into analog voltages. Analog voltages to be 
converted to digital numbers must lie in the zero to minus ten volt range. 
Similarly, digital numbers converted to analog voltages result in voltages 
between the same limits. This restriction of voltages to and from the D/A 
conversion equipment will usually require the application of bias voltages 
to the actual analog system quantities. 

For analog computer simulation studies, the necessary biasing is ac- 
complished by the application of a five volt bias to the actual system 
voltages in an operational amplifier used as a summer. The inherent sign 
inversion will result in the proper input to the ADC. For example, if X 
represents the analog quantity to be converted to a digital number, then 
the actual input to the.ADC is -(X +5) which will lie in the conversion 
range if [|X| < 5 volts. Table 2-1 is a analog voltage-digital number 
conversion table for the ADC-DAC package. It should be noted that positive, 
unbiased voltages result in positive digital numbers (bit 11 is "O"). The 


same sign correspondence is also applicable for negative values. 


Unbiased Biased Digital 
Analog Analog Number 
Voltage Voltage (octal) 
(volts) (volts) 


4000 
4632 
5463 
6315 
7146 
0000 or 7777 
0631 
1462 
2314 
3144 
3777 





Table 2-1. Analog Voltage-Digital Number Conversion Table. 
Approximately 100 microseconds is required for each analog-to-digital 


number conversion and approximately 20 microseconds.’ is required for digital- 





to-analog conversion. 

Accessory equipment included with the DAC-ADC package includes a 
common five-volt bias source, a "input disable" jack and a relay which is 
controllable by means of an external function code from the CDC 160 computer. 
The relay operation is explained in Appendix II. The common bias source pro- 
vides a common bias to all ADC inputs and DAC outputs thus eliminating in- 
div iaena bias adjustment problems and errors. 

The..““inputidisable" jack is connected by means of an AND-gate to the 
"input ready" line of the CDC 160 computer input capable and provides a 
means of delaying the CDC 160 computer by an external timing device. To in- 
put a digital word from the ADC, the CDC 160 computer sends an "input request" 
to the ADC at which time the ADC converts the analog voltage to a digital 
number on the selected A/D channel. Upon completion of the conversion pro- 
cess, the ADC sends an “input ready" signal to the computer after which the 
computer will input the digital words. If, however, the “input ready" is 
held at ground level by means of an external device, the “input ready" 
signal to the computer is delayed until the external device drops the "in- 
put disable" to -13 volts. Thus, external sample timing control may be 
accomplished via the “input disable" jack. Note, however, that when the 
"input ready" line is released by the external clock, a stale analog sample 


value will be transmitted to the digital computer and should be discarded. 


2.4 PACE TR-20 ANALOG COMPUTER. 

The Pace TR-20 analog computer is a solid state analog computer employ- 
ing twenty operational amplifiers. Saturation levels on these amplifiers are 
(+) ten volts. Standard analog computational techniques are applicable to 
this computer. External connection to the “operation-reset" relays is pro- 
vided which permits remote control of the computer from the CDC 160 digital 


computer using the accessory relay of the ADC-DAC equipment (See Appendix II). 
7 
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CHAPTER 3 
GENERAL HYBRID CONTROL SYSTEM CONSIDERATIONS 
3-0. GENERAL LIMITATIONS IMPOSED BY HYBRIDIZATION. 
If a digital computer is to be used in the control of an analog system, 
then the following factors may, in general, affect the system response: 

A. The inherent time delay required to convert a sequence of 
analog quantities to digital numbers. 

B. The inherent time delay imposed by the serial nature of compu- 
tation in a digital computer (as opposed to parallel computation in an ana- 
log sense). | 

C. The manner in which computational constants are entered into 
the computer. z 

D. The method used to obtain sample timing control if the analog 
inputs to the digital computer are to be sampled at regular intervals. 

E. The inherent magnitude limitations imposed by the A/D con- 
version process. 


Each of these factors must be considered in the implementation of hybridiza- 


tion, and are discussed more fully below. 


3.1 A/D CONVERSION TIME DELAY. 

The analog-to-digital conversion of an analog voltage cannot be done 
instantaneously but requires a finite time dependent upon the type of con- 
version process used and the specific equipment package. As explained in 
Chapter 2, the successive approximation type converter used for this project 
requires about 100 microsensiftisi for the conversion process. In many appli- 
cations, a simultaneous sampling of several analog quantities is desired; 
however, due to this finite conversion time, simultaneous sampling cannot 
be achieved by the A/D converter. The seriousness of this problem is deter- 


mined by the number of "simultaneous" samples desired and the accuracy of 
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solution desired. 

As will be explained in Chapter 4, four channels of A/D input were 
sampled in succession for a particular simulation problem with negligible 
effect noticeable due to non-simultaneous sampling. However, if a large 
number of inputs are desired, time delay effects may become significant. 

In this case, an external samp le-and-hold device will undoubtedly be neces- 
sary to sample all inputs simultaneously and hold the individual samples 
until the ADC completes conversion on all channels. 

3.2 COMPUTATION TIME DELAY. 

A digital computer aie necessarily perform arithmetic operations in a 
step-by-step or serial manner. When used for hybrid control purposes, the 
usual application involves a series of arithmetic operations on sampled in- 
put quantities and the “outputting” of the result. Due to the serial nature 
of computation, a finite time exists between input of samples and output of 
result. This time delay is a function of the cycle time of the specific 
computer and the complexity of computation. 

Efficient computer programming is obviously required to minimize these 
effects and in some applications a prediction scheme may be called for if 


extreme accuracy is desired. 


3.3 ENTRY OF CONSTANTS INTO THE COMPUTER. 

Regardless of the type of computation required in a hybrid control 
computer, certain constants must be entered into the digital computer. If 
the computer has a floating point capability, constant entry does not pose 
a serious problem. However, for a fixed point machine (such as the CDC 160), 
the format in which these constants are entered into the computer may have a 


¢ 


pronounced effect on the amount of programming required. 


Appendix I explores this problem in detail for the CDC 160 computer 
which requires access to an external arithmetic unit for multiply and divide 
operations. To minimize the computation time delay (Paragraph 3.2) and to 
provide ease of programming, the method of constant entry used herein is to 
convert all constants to the number system of the analog-digital converter 
(Table 2-1). All arithmetic computations are then done in this "psuedo- 
number" system. See Example 4, Appendix I. Specifically, the method used 
herein is: 

A. Enter constants into selected storage cells of the 160 com- 


puter using normal decimal to octal conversion procedures to obtain the 


octal representation of a feciael number. The format of entry is (x ERS,» 
i.e., the radix point is implied after the most significant octal digit. 
The number (2-5)46 would be entered as (2400) ,. 

B. Conversion of the constants from normal octal notation to the 
"pysuedo-number" system of the ADC (Table 2-1). This is done by means of a 
subroutine (See Appendix III, Subroutine KMOD), and is accomplished by multi- 
plying the constant to be converted by (0631), which corresponds to unity 
in the psuedo-number system. For example, the number (2.5) 4 = (2400) , is 
converted to (1776) ,. 

C. All further arithmetic operations are done in this new number 
system. 

As shown in Example 4, Appendix I, every multiplication must be followed by 
division by unity. This seems somewhat cumbersome, but allows the use of a 
wide range of constants and a minimum of programming to obtain a single 12- 
bit number from the 168 arithmetic unit. 

The above method of constant entry is tailored specifically for the CDC 


160 computer. The method of attack on this problem may be quite different 


for a different digital computer. Figure 3-1 is a flow chart for a general 
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program which computes an output quantity which is a linear combination 
of two sample values. It should be noted that constant conversion is re- 


quired only once and may be done outside of the sample loop. 


3.4 SAMPLE TIMING CONTROL. 

The digital computer will normally be required to sample the inputs at 
regular intervals. A method must be established to obtain this sample inter- 
val and can be done either internally or externally. Internal timing can 
be programmed into the computer using the known memory cycle time of the 
computer and having the computer perform some fixed time operation repeated- 
ly until the desired time delay is accomplished. The "time delay loop” is 
entered after the desired computations are completed in each sample cycle. | 
(See Figure 3-1). This method suffers from two disadvantages. First, it is 
a rather inefficient method ae aciee of a digital computer and secondly, the 
actual time delay is dependent upon variations in computation time of. arith- 
metic operations performed on the sample inputs, i.e., if the amount of com- 
putation varies from cycle to cycle, then the actual time interval between 
successive samples will also vary. 

Alternatively, an external timing device can be used to either "hang 
up'' the computer (such as the "input disable" explained in Chapter 2) or 
allow time-sharing of the computer with other tasks. External timing will 
provide the most accurate timing reference since the sample timing is then 
independent of the computation program execution time. 

Internal timing control was used for the simulation scheme used in 
Chapter 4, however, provisions were included in the program for external 


timing control, if desired. 


3.5 MAGNITUDE LIMITATIONS IMPOSED BY THE ADC. 
Associated with each ADC is a magnitude Limitation imposed on the analog 


ll 


quantities which must undergo conversion. This is a physical equipment 
limitation and will depend upon the specific ADC-DAC conversion package. 

This limitation not only limits the magnitude of the analog quantities to 

be converted, but also limits the arithmetic computation range in the digital 
computer if all arithmetic operations are carried out in the analog-digital 


psuedo-number system. This problem is covered in detail in Chapter 4. 
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CHAPTER IV 
SIMULATION OF AN OPTIMAL DISCRETE-TIME CONTROL SYSTEM 
4.0 GENERAL 
The presence of an analog-digital interface in a control scheme may 
have a pronounced effect on tye habnie response of the system under cer- 
tain conditions. These effects are due primarily to the A/D-D/A conver- 
sion process and must be taken into account if successful control is to be 
accomplished over widely varying conditions. This chapter will describe 
the simulation of a specifit system in order to illustrate the possible 
limitations imposed by hybridization. The system to be simulated —. 


ed by the ordinary differential equation 
X = x 4 lu (4.1) 


X(0) = C¢ 
The CDC 160 computer will be used to compute the forcing function, U, based 


on the sampled values of the state variables X, and X,: The simulation 


1 
will be accomplished in terms of both the regulator problem and the non- 


autonomous problem. 


4.1 THE REGULATOR PROBLEM. 
Bertram [1] has shown that for a sampling interval of one second, a 
time-optimal forcing function, U, for the above system is 


U(k) = - X, (k) = 1. 5X, (k) (4.2) 


where X, (k) represents the value of the state variable X. at the kth sampl- 
ing instant and U(k) represents the control value to be applied to the sys- 
tem between the kth and k+lst sample. This control law will force all of 
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the state variables to zero in minimum time. The system is shown in block 
diagram form in Figure 4-1. 

The continuous plant will be simulated on the PACE TR-20 analog computer. 
The CDC 160 digital computer will be used to sample the state variables X, 
and Xo at one second intervals via the ADC and compute the forcing function 
U(k) which then will be G¥eenitted to the continuous plant via the DAC. 


The DAC acts as a zero-order hold for the continuous plant. A block dia- 


gram of the hybrid system simulation is shown in Figure 4-2. 


4.1.1 LIMITATIONS IMPOSED BY HYBRIDIZATION, 

Simulation of the analog aoatee of the system is done in a conventional 
manner, the only restriction being that the maximum analog values of the 
state variables lie within (+) five volts (prior to biasing) so as to remain 
within the voltage limitations imposed by the ADC. This limitation, at 
most, will impose a conventional analog computer magnitude scaling require- 
ment. The bias voltages must be applied to ensure that the analog state 
variable voltages lie within the zero to minus ten volt range for A/D-D/A 
conversion. 

The CDC 160 computer program appears to be relatively simple and straight- 
forward following the flow chart shown in Figure 3-1. However, if a wide 
range of initial conditions and state variable values is to be consideréd, 
then additional limitations impose themselves on the programming of the 
computer. Consider, for example, the formation of the product 


C12) guctinling ten GOs 0 


The 160 computer generates this product as 


(6632), x (3144), = (7153). = (-1) 46 


This inconsistency can be attributed to the fact that the magnitude of the 
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Figure 4-2 Block diagram of Simulation Scheme 
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desired product is greater than five while the 'pseudo-number system" 
associated with the analog-digital interface is consistent only for numbers 
of magnitude less than five. The desired product (-6) has "overflowed" the 
number system and erroneous results occur. In general, for the product 
aX, overflow occurs unless 

IX, | ee Ce (4.3) 
A similar type of overflow occurs when a sum of a series of numbers exceeds 
the overflow limit, a. 


< 
aX) + aX, soota X ==) (4.4) 


This overflow problem can be attacked in one of two ways. First, a 
suitable scale factor may be applied to the analog values such that the 
maximum scaled values of the state variables will not cause overflow. 

This method is simple but can lead to inaccuracies as the state variables 
tend toward zero where the analog-digital conversion errors and digital 
truncation effects become significant with respect to the magnitude of the 
scaled state variable values. 

Secondly, normal scaling may be used, but each product and sum can be 
tested for overflow as it is formed, and if overflow is detected, then the 
computer can be programmed to substitute the maximum permissible value (+5) 
for the intended product or sum. Quite obviously, if an optimal control 
scheme is being used, the trajectory will no longer be optimal if overflow 
occurs, but the steady state accuracy will not be affected once the state 
variables are finally driven to zero. This second approach will be develop- 
ed herein since this method must be incorporated into the CDC 160 computer 
program. 

The effects of this overflow problem on the specific system being -simulat- 


ed may best be illustrated on the X, vs X, phase plane. The area of concern 
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on this phase plane is the area bounded by the maximum digital-analog 
conversion limits, i.e., a square area centered at the origin bounded by 
|X, | = 5 and |X, | = 5. Overflow boundries for each of the products may be 
computed from Equation 4.3. Thus 


Ix, | <5 (4.5) 


and 


\X,| 5 5/25 = 3.33 (4.6) 
The summation overflow boundries may be computed from Equation 4.4, i.e., 


|-X, - 1.5X £5 


1 2| 
or e 
-.667X, - 3.33 X, S -.667X, + 3.33 (4.7) 
The six straight lines defined by the equality relationships of Equations 
4.5, 4.6 and 4.7 reduce the original digital-analog conversion area of the 
X, Vs Xo phase plane as shown in Figure 4-3. In this figure, Region I 
represents the product overflow area; Region II the summation overflow 
area; and Region III represents the area where both product and summation 
overflow occurs. Any initial condition or state variable sample value which 
lies in any of these regions will cause an overflow condition with result- 
ant errors in computation. It should be noted that these overflow bound - 
ries are determined directly from the values of the state variable feedback 
coefficients, a. If these coefficients are changed, the overflow bound- 
ries also change according to Equations 4.3 and 4.4, 
The dynamic behavior of the system described by Equation 4.1 can be ex- 
pressed in terms of the discrete difference equation 


1 1 % 


X(k+1) = X(k) + U(k) (4.8) 
0 1 1 
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Region I 


Region III 





Figure 4-3. Phase Plane Overflow Regions 
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For a time-optimal trajectory, the control input must satisfy Equation 4.2. 


Combining Equations 4.2 and 4.8, the difference equation can be reduced to 


5 
X(k+1) = X(k) (4.9) 
-1 -3 
It is important to note that Equation 4.9 is valid only if X(k) lies within 
the overflow boundries shown in Figure 4-3. If the state point lies out- 
side of the non-overflow region, U(k) will not take on the proper value for 
optimal control, hence Equation 4.9 no longer holds. 
Keeping the above Mintte clone on Equation 4.9 in mind, not only must 
the state point lie within the overflow boundries for the kth sampling 
instant but also at the k+lst sampling instant. In other words, the kth 


control input must not drive the state point into the overflow region. 


From Equation 4.5, 4.6 and 4.9. 


JX, (k4+1)| = [2X,(k) + 2X (k)| S 5 
or 

-2X,(k) - 20 s X, (Ie) s -2X,(k) + 20 (4.10) 
and 

[x5 Cee) | = |X) a | 
or 


-2X,(k) - 6.66 < X, (Ie) < -2X, (k) + 6.66 (4.11) 


The two limiting straight lines expressed in Equation 4.10 lie entirely 
outside the digital-analog conversion boundry and impose no restrictions; 
however, the two lines defined by Equation 4.11 further reduce the non- 
overflow area of the phase plane as shown in Figure 4-4. This figure shows 
the optimal control boundries of the phase plane; if the initial state 
point is located within this boundry, optimal control will be realized. 
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It should be emphasized that the state point may be driven to the 
origin from any point within the digital-analog conversion boundry provided 
that the CDC 160 computer program is written to correct for overflow; how- 
ever, the trajectory will be time-optimal only if the initial state point 


lies within the boundries shown in Figure 4-4. 


4.1.2 PROGRAMMING TO DETECT /CORRECT FOR OVERFLOW. 

The CDC 160 computer can be programmed to detect and correct overflow 
conditions. The method of detection of an overflow condition is dependent | 
upon the operation Oya iret multiplication) causing the overflow. Re- 
ferring to Figure 4-5 (a schematic representation of the digital-analog 
psudeo-number system), it is seen that the summation overflow path is con- 
tinuous across the overflow boundry and a change of sign always occurs for 
a summation overflow. Note also that overflow as a result of the addition 
of two numbers only occurs if both of the numbers agree in sign. Thus sum-— 
mation overflow detection can be programmed in the following manner. 

1. Compare the signs of the two arguments. 

2. If the arguments differ in sign, overflow cannot occur. 

3. If the arguments agree is sign, overflow is possible. In this 
case, if the sign of the resultant sum acreectethtthe Bign Cle 
arguments, overflow has not occurred. If the sign of the sum 
differs, overflow has occurred. 

The product overflow path, on the other hand, cannot cross the overflow 
boundary as can be seen in Figure 4-5. The overflow product does not change 
sign but undergoes a discontinuous jump in magnitude as overflow occurs. 
Product overflow can only occur if both of the arguments are greater than 
unity. Thus to detect overflow of the product aX: 

1. Determine if a. is less than unity. If so, overflow cannot 


occur and no further tests are required. 
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2 If as is greater than unity, compare the magnitudes of 


aX, and Xo If the magnitude of the product is less than the 
magnitude of Xx. itself, overflow has occurred; otherwise not. 

Overflow correction in both cases consists of substituting the maximum 
possible value (+ 5) for the overflow quantity consistent with the signs of 
the arguments causing overflow. 

To detect and correct for overflow in the system being simulated, a 
product overflow test must be performed on the product (-1.5X,) and a sum- 
mation overflow test performed on the sum (-X,) ~~ (-1.5X,). If the program 
is written to allow variacieh of the feedback coefficients, then the pro- 


duct a,X, must also be tested to allow values of a, greater than unity. 


i 1 1 


-_ 


4.2.0 THE NON-AUTONOMOUS PROBLEM. 

Suppose it is desired to force the system described by Equation 4.1 
to follow an external signal of the form 

R(t) = r,t rt + (4.12) 

If one associates X, (t) with system position and X(t) with system velocity, 
then R(t) can be considered as "target" position and r as "target velocity. 
Bertram's time-optimal solution of the regulator problem (Equation 4.2) is 
designed to drive the state variables of the original system to zero in 
minimum time. This same optimal solution can be applied to the non-auto- 
onomous case if one defines a position error and a velocity error and con- 
siders these errors as the state variables of a new system. With the inputs 
thus "imbedded", the time-optimal solution will drive these errors to zero 
in minimum time. 


Define: 


i 
i 


position error = ¥, (t) X, (t) -R(t) 


I 
i 


velocity error ¥,(t) X(t) = (4.13) 
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Applying Bertram's time-optimal forcing function, U, to this new system: 


U(k) = -Y, (k) - 1.5Y,(k) 
which by substitution of Equations 4.12 and 4.13 becomes 
U(k) = -[X,(k) - r, - ft(k)] - 1.5[X,(k) - £] (4.14) 


The simulation of this system can be accomplished in the same manner as 
the regulator problem adding the requirement that the CDC 160 computer 
take into account target position and velocity in addition to system 
position and velocity. The simulation scheme chosen is shown in Figure 
4-6. The value r, is set into the CDC 160 computer manually due to the 
four-channel Heteseiecon the ADC. Note that in an actual system, it is 
possible for the signal R(t) to be transmitted to the digital computer 


from an entirely remote source such as a radar system. 
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4.2.1 LIMITATIONS IMPOSED BY HYBRIDIZATION, 
As with the regulator problem, the most severe limitations imposed 
on the non-autonomous system are due to overflow phenomena which impose 
maximum limits on the signal R if optimum control is to be realized. In 
the computation of Equation 4.14, the 160 computer must perform four addi- 
tions and two multiplications. At first glance, it appears that a total 
of six overflow tests must be performed; one for each operation involved 
in the computation of U(k). Closer examination will show that the actual 
number of overflow tests required is the same as for the regulator problem. 
Equation 4.14 can be written as: 


_ Uk) = a5, (k) + a5, (k) (4.15) 


where 


Sy) a eres, 


S, (k) = X, (k) - fr 


a, = -l for optimal control 


-1.5 for optimal control. 


fo 
ll 


The 160 computer must form S, in two steps. If the first step is the 
addition of (-r) to KX, (k), then this sum may possibly overflow. If over- 
flow does occur, reference to Figure 4-5 will show that summation overflow 
is continuous across the overflow boundry. Further addition of the eeencies 
(-rt(k)), if of appropriate sign and magnitude, may thus be able to bring 
the entire sum Ss) out of the overflow condition. If it cannot, then S 


itself overflows. The same result obtains if the order of addition is 


1 


reversed. Therefore an overflow test need not be performed on the inter- 


mediate sum in the computation of S, and only one overflow test is required 


1 


on the entire composite sum. 
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Noting that for a one second sampling interval, t(k) = k, then Equa- 
tions 4.14 and 4.8 may be combined to yield: 
a 45 


X(k+1) = X(k) + [r + i{k + roy) (4.16) 
-]) S2|~ 1 


Assume X(0) = 0. Then from Equation 4.16, 


5r + .75r 
Oo 


X(1) = (4.17) 
r + 1.5r 
Oo 
. ie kr 
X(k) = + 5 Kos 2 eo, ; (4.18) 
- 
To prevent overflow of the sum 54: 
= | ; - P 
js, Qe)| = |X, (k) etd || SS) = (4.19) 


Evaluating this inequality for various values of k by the substitution 


of Equations 4.17 and 4.18 (still under the assumption that X(0) = 0), one 


obtains: 
CEs |e= seterek = 6 (4.20) 
ee Oc r << 20- ar, for k = 1. (4.21) 


(3) Is, (KD | = 0 fork 2 2. 


Equations 4.20 and 4.21 define the overflow boundries of s5 as a function 


of rand r. S. can only overflow for k = 0 ork = 1. (Note that these 


1 
same limitations on r. and r also hold for the quantity a,S,(k) if a, = 


1 
-1). 
Similar overflow analysis can be performed on the quantities S,(k); 
aS, (k) and U(k). Table 4-1 lists the results of this overflow analysis. 
If the straight lines defined by the equality relationships of each of the 


input limitations of Table 4-1 are plotted on the rvs r plane, the over- 


flow boundries for the system are defined. The most severe limitations are 
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determined by the quantities a and U(k), and thus determine the range 


2°2 
of inputs allowable for optimal control. These limiting boundries are 
shown by the solid lines of Figure 4-7. [It should be noted that Table 
4-1 and Figure 4-7 are valid only wonder the assumption that X(0) = 0. 
For this special case, the optimal control region under non-zero input 
conditions is equivalent to the optimal control region on the x) vs Xo 
phase plane for the regulator problem. (See Figure 4-4). 








INPUT LIMITATIONS 


S 3.33 - -667r_ 


ae WE ae 6674 





Table 4-1. Input Limitations for X(0) = 0. 
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A further limitation on problem running time is imposed due to the 
fact that X) will eventually exceed the five volt A/D conversion limit as it 


follows target position. This limitation can be expressed as: 

je["S 5-2 XCOVFMPO, & = 1a. (4.22) 

, i 

The limiting case of this relation for various values of k is shown by the 
dotted lines of Figure 4-7. If a problem running time of ten seconds 
(k = 10) is desired, for example, then|r | < 5 r/10. In a practical 
situation, A/D conversion equipment must be chosen so as to provide a 
conversion limit which exceeds the maximum expected value to be converted, 
or alternatively, successive scale factors must be applied to the quantity 
X) as the conversion limit is approached. This limitation cannot be remov- 
ed by appropriate programming of the digital computer. 

As a result of the above overflow analysis, and if one programs the 
digital computer to permit variation of the feedback coefficients ay and 
ov) (and thereby varying the overflow boundries) ,the following quantities 
must be tested and corrected for overflow: 

(1) The product a,5, (k) =a) [X, (k) Be rt(k: ]. 

(2) The product a8, (Kk) = ao [X,(k) - rj. 

(3) The sum U(k) = a,5S, (k) + aS. (k). 
If ty and r are zero, then the above equations reduce to the overflow 
quantities considered in the regulator problem discussed in Section 4.1. ° 
Thus one computer program will permit control of the system under both zero 
and non-zero input conditions. Program DIGIT (Appendix IV) is a digital 
computer program written for the CDC 160 computer which performs the above 


overflow tests and corrects the overflow condition based on the simulation 


scheme shown in Figure 4-6. 
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(1) U(O) overflow boundry . 
(2) -1.5S,, (0) overflow boundry. 
(3) -1.58',(1) overflow boudry. 


--- X, overflow at time indicated. 


Figure 4-7. Input Limitations for X(0) = 0. 
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4.3 SIMULATION RESULTS, 

The system described by Equation 4.1 was simulated on the CDC 160 
digital computer and the PACE TR-20 analog computer using the simulation 
scheme shown in Figure 4-6. Figures 4-8 thru 4-13 are actual x-y recorder 
plots of the system response for various initial conditions and input 
signals. Figure 4-10 shows system response under overflow conditions. It 
is seen that two eee ae required to drive the state point into the 
optimal control region after which optimal control is exercised. For all 
plots, the program DIGIT internal timing loop was used. The finite response 


time indicated on the signal U(t) is due to recorder lag. 


“" 
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Simulation Results for Initial Position Error 


| —__- | > * ©) 
3 





X(0) 
5 


=a 
{ 
i 


ii 
© 


R(t) 


>t 








} 
er fp ep ef et 
1 2 3 4 
“14 
~2 2. 


Figure 4-9 Simulation Results for Initial Position and Velocity Error 








Simulation Results for Initial Position Error 
Resulting in Non-Optimal Control. 
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Simulation Results for Step Input 
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Figure 4-12. Simulation Results for Ramp Input 


39 





Xo(t) x(0) = 0 
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Simulation Results for Step plus Ramp Input, 
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CHAPTER 5 
SUMMARY 

One of the first problems encountered in the use of a short word 
length digital computer for hybrid control purposes is the problem of 
entry of computational constants into the computer in a format compatible 
with the A/D conversion process. Two choices present themselves at this 
point; namely, (1) to perform all arithmetic operations in the A/D psuedo- 
number system or (2) to perform all arithmetic operations in the standard 
binary number system normally used by the computer. In the first case, 
the computational sctesaat dude be converted to the A/D number system and 
in the second case, the digital sample values must be converted to normal 
binary notation. The first choice was chosen for this project resulting in 
subroutine KMOD which can be used for a large class of hybrid control pro-. 
grams wherein constants must be entered into the computer manually and must 
be capable of being changed by the program user. 

The above choice of number system results in the phenomena of number 
system “overflow'’ within the computer during the computation of the control 
law. As demonstrated in Chapter 4, each term of the control law is capable 
of overflow under certain conditions. The two types of overflow (product 
Or summation) are characterized by certain properties which allow detection 
of an overflow condition by the digital computer, and if overflow is de- 
tected, the computer can be programmed to partially correct for this over- 
flow condition. This overflow detection/correction programming can greatly 
increase the range of initial conditions or system inputs even though optimal 
control is not realized under overflow conditions. 

Program DIGIT, written for the specific system simulated in Chapter 4, 


is sufficiently general so that it may be used with any second-order analog 
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simulation where the control law to be computed by the digital computer 
consists of a linear combination of the state variable values. A wide 
range of feedback coefficients and variable sample timing can be accommo- 
dated by this program. 

Further investigation is needed on the second method of arithmetic 
computation, i.e., the conversion of the A/D digital samples to normal 
binary notation. It is believed that this method is most easily imple- 
mented on a long word length computer possessing the capability of float- 
ing point arithmetic. The programming difficulties presented by this 
conversion process appear 8 be excessive for a short word length, fixed 
point arithmetic machine such as the CDC 160 computer, and since this 
conversion must be included in the sample loop of the program, excessive 


computation time between the input of samples and the output of the con- 


trol value may result. 
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APPENDIX I 


OPERATION OF THE CDC 168 ARITHMETIC UNIT 
I-1. MULTIPLICATION AND DIVISION OPERATIONS. 

The CDC 168 arithmetic unit is used in conjunction with the CDC 160 
computer for single precision (short) multiplication and division in the 
eae developed herein. A detailed discussion of these two modes of 
operation is given below. 

A. Single Precision (Short) Multiply. To initiate the short 
multiply operation, the 160 computer selects the arithmetic unit and trans- 
mits a 12-bit multiplicand atl a 12-bit multiplier to the arithmetic unit. 
The steps by which the arithmetic unit computes the product are shown in 
Figure I-1. Upon completion-of the short multiply operation, the product 
is located in two successive storage locations in the 160 computer; the 
least significant half (LSH) in the first storage cell, the most signifi- 
cant half (MSH) in the second cell. 

B. Single Precision (Short) Divide. To initiate the short 
divide operation, the CDC 160 computer selects the arithmetic unit and 
transmits two 12-bit words which form the dividend followed by a 12-bit 
divisor. The two 12-bit words forming the dividend must agree in sign. 
Figure I-2 shows the step by step operation of the arithmetic unit in the 
short divide mode. It is to be noted that the 160 computer transmits the 
least significant half of the two word dividend first, and inputs the 12- 
bit quotient followed by the 12-bit remainder. 

I-2. SHORT MULTIPLICATION AND DIVISION EXAMPLES, 

As noted above, the short multiply operation yields a two word pro- 

duct in the 160 computer. Successive operations in the 160 computer, and 


digital-to-analog conversion are single word (12-bit) operations. Thus a 
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From CDC 160 





Input word 
#1 and #2. 
(12 bits each) 













Remove sign 
Bit (bit 11). 
From word #1 


Sign bit #1 
= 1 or 0? 


Remove sign 
Bit (bit 11) 
From word #2. 


Sign Bit #2 =] 
= 1 or 0? 

Form 22-bit 
Product 

Sign Bit #1 No 

= sign bit #2? 











Complement 


Remaining 
bits’ of word #1 













Complement 
Remaining 11 
Bits of word #2 


Product sign 
Bits = 1 











Product sign Complement 
Bits = 0 22-bit 
product 


Form LSH of 
product 
See Note l. 





Form MSH of 


product. See 
Note l. 


Output LSH 


of product to 
160. See note 2 


Output MSH 
of production 
160. See note 2 


To CDC 160 


NOTE 1. The least significant 
half (LSH) of the final product 
is formed from the least sign- 
ificant 11 bits of the 22-bit 
product and a sign bit. The 
sign bit is inserted as the 
most significant bit, thus 
forming a 12-bit LSH. The 
most significant 11 bits of 
the 22-bit product are 
similarly treated to form the 
MSH of the final product. 


NOTE 2. Output is initiated 
by the 160 computer. 


Figure I-1. CDC 168 Arithmetic Unit Short Multiply Flow Chart 
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From CDC 160 


Input LSH of 
dividend 
(word #1) 












Input MSH of 
dividend 
(word #2). 


Input 
Divisor 
(word #3) 









Remove sign 
bits of 
words 1,2,3. 


Sign Bit #1 = 1 Complement 
= 1 or 0? 22-bit 
= 0 dividend 







Sign bit #3 Complement 
= 1 or 0? ll-bit 
Divisor 






Form ll-bit 
Quotient and 
ll-bit remainder 


Sign bit #1 No 
= sign bit #3? 


Sign bit of 
Quotient and 
Remainder = 0 













Sign bit of 
Quotient and 
Remainder = l. 











Complement 
Quotient ard 
Remainder 








Output ll-bit 
Quotient plus 


sign bit 





Output ll-bit 
Remainder 
Plus sign bit. 






To CDC 160 


NOTE 1: LSH = least signi- 
ficant half, MSH = most 
significant half. 


NOTE 2: Output execution 
is initiated by the 160 
computer 


Figure I-2. CDC 168 Arithmetic Unit Short Divide Flow Chart. 
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12-bit, single word product must be retrieved from the two word answer 


received from the 168 arithmetic unit. The manner in which this reduc- 
tion in word length is accomplished is determined by the manner in which 
the arguments of the multiply operation are represented in the 160 computer. 
Assume that one of the arguments (either multiplier or multiplicand) will 
result from the analog-to-digital conversion of a voltage, and hence will 
be of the format shown in Table 2-1. Assume also, that the second argument 
will be a constant entered into the computer by the user. Then a choice 
exists as to the manner in which this constant is represented in the 160 


computer. Let 


and 


ma = )0. 1, 2,:45., © Ob. 
If the symbol K represents the constant to be entered, then K can be 
entered in one of the following ways: 

1. As an integer, i.e., 


K = (nn...),, = (mmm ) 


10 
2. As a fraction, i.e., 
K = (O.nnn), > = (.mmmm) . 


3. As a combination of the above, i.e., 


K 


(n.nnnn), 5 = (m.mmm) . 
or 


K 


(nn.nnn), = (mam .mm) ., 
The radix point is not physically present in the computer word, but is 
included above for clarity. 

A. Example 1. Consider the product, (K) x (Y), where K = (-2) 46 


and Y = (+2) 56 = (1462), after analog-to-digital conversion. Using method 
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(1) to enter K into the 160 computer, then (+2) 16 = (0002) , and to 


represent (-2)50> the 7's complement is used, i-e., 
K = (-2) 56 = (0002) . = (7775). 


Referring to Figure I-l, the 168 arithmetic unit will form the product (K) 
x (Y) as follows: 


[(7775), - sign bit] = (11 111 111 101), (sign bit = 1) 


(7775), - sign bit|= (00 000 000 010), 

[(1462), - sign bit] = (01 100 110 010), (sign bit = 0) 
(22-bit product) = (0 000 000 000 011 001 100 100), 
(22-bit product) = (1 111 111-111 100 110 011 011), 
Addition of the sign bits to form the LSH and MSH yields: 


(4633), a, 


LSH of (K) x (Y) (100 110 011 001), 


MSH of (K) x (Y) (111 111 111 111), (1777). (-9),5 6 


Thus, if the constant K is entered as an integer, with the radix point 
(implied) located to the right of the least significant octal digit, then 
the desired 12-bit product is located in the LSH of the two word product 
obtained from the arithmetic unit and no further operations are required to 
retrieve the 12-bit product. 

B. Example 2. Consider the product (K) x (Y), where K = (0.25), 


and Y = (+4) = (3144), after analog-to-digital conversion. Using method 


10 


(2) above to enter K into the 160 computer: 


K = (0.25),, = (2000). 


The 168 arithmetic unit will form the product (2000) , x (3144) , as follows: 


0) 


li 


(2000), - sign bit = (01 000 000 000), (sign bit 


(3144), - sign bit = (11 001 100 100), (sign bit = 0) 
22-bit product = (0 110 011 001 000 000 000 000), 
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Addition of the sign bits to form the LSH and the MSH yields: 

LSH of (K) x (Y) = (000 000 000 000) = (0000). = (0), 9 

MSH of (K) x (Y) = (001 100 110 010), = (1462), = (+2)) 6 
If the MSH is now right- shifted one binary bit, i.e., 

MSE Wars ( Es) eye — 00081 Omari 001), = (+1)) 9; 
then we have the desired 12-bit answer in the MSH of the two word product 
received from the arithmetic unit. Further difficulties arise using this 
method, however, due to the use of bit 11 as the sign bit in the 160 com- 


puter word. For example, *if K = (+0. 5), then the octal representation 


of Kis (4000).,. Both the computer and = arithmetic unit interpret this 
number as a negative number since bit 1l is a''l''. Hence the usable posi- 
tive fractions are restricted to be less than one-half. 

C. Example 3. Consider the product (K) x (Y), where K = (2. 16 
and Y = (2) 16 = (1462) Using method (3) above to enter K into the 160 


computer, 


eee 5)16 = (2400)... 
The arithmetic unit will form the product (2400), x (1462), as follows: 


0) 


4)- 


(2400), ~ sign bit = (10 100 000 000), (sign bit 


(1462), - sign bit = (01 100 110 010), (sign bit 
22-bit product = (0 011 111 111 101 000 000 000), 


0) 


Addition of the sign bits forming the LSH and the MSH: 
LSH of (K) x (Y) = (001 000 000 000), = (1000), 


Msi of (K) =x (Y) —(000 aie 111), = (0777), 


At the completion of the multiply operation, the LSH is located in storpge 


cell n and the MSH is located in cell n+l of the 160 computer, i.e., 


ay 





FC Del i Pm a lm Nae i a 
.000| 111] 111] 111 |f 001 | 000 | 000 | 000° 
Cell n+l Cell n 


The desired answer (3775). can be uncovered" by the removal of bit 11 of 
the LSH (the sign bit added by the arithmetic unit) and a "long right inter- 


cell" shift of one binary bit in the two cells: 
000 111 111 111 001 000 000 000 
vee 4 ™“™~, 


YN 
000 | 011 [111 |111 |f 101 | 000 | 000{ 000 
0 fee eT | 5 Ooms 
(MSH) (LSH) 


However, the answer is split between two storage cells. A series of shift- 
ing and masking operations must be performed on these two words to place 
the desired answer in one 12-bit storage cell. This problem could be 
greatly simplified by the incorporation of a long right shift capability 
in the 168 arithmetic unit, however, at present, this shifting must be 
done in the 160 computer. 

D. Example 4. Consider the product (K) x (Y), where both K and Y 
are expressed in the number system shown in Table 2-l, i.e., in the analog- 
digital conversion number system. Let K = Y = (+2) 56 = (1462) ,. Then the 


arithmetic unit will form the product (K) x (Y) as follows: 


LSH of (K) x (¥) = (010 111 000 100), = (2704), 


MSH of (K) x (XY) = (000 101 000 110), = (0506), 


The desired answer (3144), does not appear in either the LSH or the MSH. 
However, if this two word product is divided by (0631), (which corresponds 

to unity in Table 2-1), then the desired product will appear as the quotient. 
Referring to Figure I-2, let the dividend by (05062704) , and the divisor be 
(0631)... Then the arithmetic unit performs the short division in the follow- 
ing steps: 


50 


~ 





ahem 


veh & 


i. 

















? 





a —_ 





és 


eave al 
" aaa 7 a) ee oa Ce 8) 21 eee 
we ct dm gl —- - 














. 


oo 





CS St ——« oe 


















(2704), - sign bit = (10 111 000 100) (sign bit = 0) 
Zz 
dividend 
(0506) , - sign bit = (00 101 000 110), (sign bit = 0) 
divisor — (0631), - sign bit = (00 110 011 001), (sign bit = 0) 


11-bit quotient = (11 001 100 100), 
11-bit remainder = (00 000 000 000). 


Addition of sign bits to quotient and remainder yields: 


quotient = (011 001 100 100), = (3144), = (+4), 


remainder = (000 000 000 000), = (0000) , = (9) 56 


By means of as short multiply operation followed by a short divide opera- 
tion, the product (2) 10 x (2) 46 has been accomplished with the desired 
answer in a single 12-bit storage cell in the 160 computer. Note that this 


computation can be expressed as: 


(1462) x (1462) = (1462) x [(2) x (0631)] 
(0631) (0631) 


Thus if all arithmetic operations involving quantities which have 
undergone A/D conversion are done in the "analog-to-digital psuedo-number 
system" as shown in Table 2-1, then the CDC 160 programming procedures are 
greatly simplified. Any constants entered into the computer by the user 
must either be entered into the computer in the psuedo-number system for- 
mat (Table 2-1), or must be converted to this format by appropriate pro- 


gramming. This latter method was chosen for work done in this paper. 
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APPENDIX II 
EXTERNAL FUNCTION CODES FOR THE CDC 160 COMPUTER 
II-1. EXTERNAL FUNCTION CODES. 
The CDC 160 computer machine language code for selection of external 


equipment is: 


CELL , CONTENTS OSAP MNEMONIC 
n 7500 EXF 
n+l XXXX XXXX 


+ 


where the 12-bit operand, xxxx, is known as an external function code. 

The external function command as shown above is used to select an external 
device to perform some specific function depending upon the external func- 
tion code used. A partial list of external function codes is given in 
Table II-1l. 

Only one external device may be selected at any one time. Selection 
of any device automatically disconnects any other selected device. If an 
illegal selection is attempted the computer will be indefinitely delayed. 
Appropriate input-output instructions, if required, must follow the external 


function command. 





External Function External Remarks 
Device Desired Function 


Code 















Short divide 
Short multiply 
Long divide 
Long multiply 
A/D conversion 
D/A conversion 


168 arithmetic unit 
168 arithmetic unit 
168 arithmetic unit 
168 arithmetic unit 
ADC -DAC 
ADC -DAC 







See note l. 
See notes 1,2. 





Table II-1. CDC 160 External Function Codes. 
NOTE 1. For the ADC-DAC external function codes, y is the channel number. 
NOTE 2. The octal number x controls a relay in the ADC-DAC which is avail- 
able for control of external equipment. As used herein, this relay is 


connected in series with the "'OPERATE-RESET" switch of the Pace TR-20 
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analog computer (as shown in Figure II-1) to permit the CDC 160 digital 
computer to control the analog computer. If the digit, x, in the external 
function code 24xy is "1", then the relay switches the analog computer to 
the OPERATE mode; if x is "0", the analog computer is switched to the RESET 
mode. The MASTER CLEAR key on the 160 computer console will also reset 


the analog computer. 





ADC-DAC Panel 


Pace TR-20 Front Panel 


Figure II-1. Connections for Remote Control of Analog Computer 
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APPENDIX III 

SUBROUTINE KMOD 
III-1. GENERAL DESCRIPTION. 

Subroutine KMOD is designed to convert a series of octal numbers to 

a format compatible with the CDC 160 computer, the CDC 168 arithmetic unit 
and the D/A-A/D convérsion equipment. Conversion is accomplished by multi- 
plication of the number by (0631), using the CDC 168 arithmetic unit. After 
conversion, the numbers are of the same format as the output quantities of 
the ADC. Further arithmetic operations involving these numbers and A/D 
output quantities can —— be accomplished with greatly simplified program- 
ming procedures. An abbreviated conversion table is shown in Table III-1. 
Comparison of Table III-1 and Table 1 will show that KMOD acts as a "analog- 
to-digital converter'' for these numbers (designated herein as K'). As 


with analog voltages to be converted in the ADC, the numbers to be converted 


by KMOD are restricted to the range [K"| m5. 


III-2. USE OF SUBROUTINE KMOD. 
A. Calling sequence. This subroutine may be called from the main 


program with the following two lines of coding: 


Machine Code Mnemonic Remarks 
0101 SHA 01 (or PTA) Shift contents of P-register 
to the A-register 
7000 JPI 00 Indirect jump to address 


located in cell 0000. 


The subroutine generates and stores its own linkage address to the main 
program. 

B. Storage Requirements. KMOD requires 161 (octal) storage cells. 
KMOD may be loaded at any location in the 160 memory. 

C. Data Entry. The following data must be entered into the CDC 160 


computer manually in the locations shown: 
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Cell Number Contents 


0000 Load address of subroutine KMOD. 


0001 Number of constants (K') to be converted. 
» (24 maximum). 


0002 Storage address of first converted 
constant, K. See note l. 

0003 | Sign of K', = {0000 if K', is positive 

0001 if Ky is negative. 

0004 Kj. See note 2. 

0005 Sign of K",. 

0006 : ye 
5 2 

0007 : Sign of K',. 

etc. 


Note 1. If more than one K' is ~ be converted, then the converted 
numbers (K, ) are stored in every other storage cell beginning with the 
storage location entered in cell 0002. This permits greater ease in input- 
output operations in the main program. 
Note 2. The constants to be converted must be entered into the 160 com- 
puter in the format (X. xXx) g. Example: If K' = (1.4325), this would be 
entered into the 160 computer as (1433) .- 

Manually entered data is not destroyed during execution of KMOD. Thus, 
execution can be repeated without re-entry of data. 

D. Low core cell usage. In addition to the low core cells used for 
manual entry of data, cells 0064 thru 0077 are used during the execution of 
this subroutine. After execution of KMOD, cells 64 thru 77 are available 


for other purposes. 
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Octal Number | Octal Number 
to be converted after conversion 
K' K 
4002 
4317 
4633 
5150 
5464 
6001 
6315 
6632 
7146 
7463 

or 0000 


0314 


0631 


1145 
1462 
1776 
2313 
2627 
3144 
3460 


3775 





Table III-1. Conversion Table for Subroutine KMOD. 
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IIi-3. DETAILED DESCRIPTION OF KMOD. 

If the arguments of the short multiply operation are of the format 
given in the preceeding instructions, then, as shown in Example 3, Appen- 
dix I, the desired 12-bit product is split between two successive storage 
locations in the 160 computer. Let the Me, product as received from 
the arithmetic unit be expressed as follows: 


(binary) 
(octal) 





where x 


either 17 digit (1 or 0) 


f 
lith bit of the 22-bit product 


<< 
fl 


w and z = sign bits inserted into the 22-bit product by the 
arithmetic unit. 


The following procedures are necessary to retrieve the desired 12- 
bit product, K = (k,k ak3k,) 9° 
A. Removal of the sign nie in the LSH. Since w has been inserted in 
the middle of the desired product, it must be removed. To do 
this: 
1. Set w equal to y. 
2. Shift MSH one bit to the right (this requires an eleven 


bit left shift in the computer). 


Upon completion of step 2, the two words are: 





MSH LSH 


B. Round-off correction. The radix point of the two-word product 
is implied after octal digit Ky. If the digits (nnn) , are greater 
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el bee - 


than (444). then k, must be increased by one unit. However, 
ie ky, = (7). then this round-off correction will propagate to 


ky 5 requiring that k, be increased by one unit. If (nnn), is 


3 
less than (444) 5 no round-off correction is necessary. Thus 
the steps in the round-off correction are: 

1. Determine if (nnn), is greater than (444) .. If not, 
skip the following steps. 

2. If round-off correction is necessary, determine if kK, 
is equal A (7), - If so, add one to k.. If not, go to 
next step. 

3. Add one to ky). 

C. Removal of extraneous bits. Digits m and n have now served their 
purpose and can be discarded. To do this: 

1. Mask off m and shift MSH left three bits. 

2. Mask off (nnn) and shift LSH left three bits. 


Upon completion of this step, the two words are: 


BRIM. ..: 
MSH LSH 
3. Add LSH to MSH to get desired number K = (k, k, k, Kk): 
III-4. PROGRAM LISTING FOR SUBROUTINE KMOD. 
The following pages give the OSAP program listing for subroutine KMOD. 


Sufficient comments are inserted to make the program self-explanatory. 
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Cell Mach Symbolic 
Code 
) 0602 adn 
1 4072 atd 
2 2000 ldd 
4 0620 adn 
4 4071 std 
5 hoz ldn 
6 4064 std 
ie 0404 ldn 
10 4066 std 
11 2001 ldd 
12 4070 std 
13 2002 ldd 
14 4067 std 
bp: 2200 lat 
16 0631 
17 4074 std 
20 2164 ldi 
21 4073 std 
22 2166 ldi 
23 6206 pjf 
24 2200 laf 
25 2000 
26 4065 std 
27 2166 ldi 
30 3465 sbd 
31 4o75> h std 
32 7500 exf 
29 5201 
34 7306 out 
35 0076 
36 7205 inp 
Ball 0100 
40 6105 nzf 
4) 6004 zjf 
40 0074 ki 
43, 0076 kia 
44 4000 con 


SUBROUTINE KMOD 


‘stoadd 


°631 


one 
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Remarks 


Form link address to main program. 


Store link address. 

Load the address of KMOD. 

Form recycle address for multiple K's. 
Store recycle address. 

Load initial sign address. 

Store initial sign address. 

Load initial K' address. 

Store initial K' address. 

Load number of (K')'s to be converted. 
Store to form index. | 

Load storage address for lst converted K'. 
Store above address. 

Load 0631. 


Store as argument for short multiply. 


The preceeding instructions perform 
several book-keeping tasks prior to the 
execution of the subroutine. All of 
the manually entered data is relocated 
so that upon execution of the following 
instructions, this manually entered 
data is not destroyed. The cell numbering 
is given relative to the initial load 
address. 

Load sign of K' to be converted. 

Store. 

Load K' to be converted. 

If bit 11 is 0, jump. 

If bit 11 is 1, load the number 2000. 


Temporarily store 2000. 

Reload K! e 

Subtract 2000 from K'. See *. 
Store in output location. 
Select a.u., short multiply. 


Output 0631 and K'. 


Input two word product, LSH in cell 76, 
MSH in cell 77. 
Unconditional jump around table. 


Input-output table. 


The sign of the K' under conversion has 
been determined and stored. The two 
word product has been received from the 
arithmetic unit. 


+. Res Oye eee, 
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Cell 





Mach 


Code 


Symbolic 


i 


SUBROUTINE KMOD (con't) 


abd 
laf 
lpd 


std 
lat 
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Remarks 


* Notes If kK! >3777, then bit ll =1, 
but the sign function has been taken away 
from bit 11 since the sign of K' is entered 
as a separate word. Since the maximum 
value of K' is 5000, then subtraction of 
2000 will yield a positive number. A 
later operation adds back this qua Fis 
Load MSH of two-word product. 

Mask off bits O01 thru ll. 

If bit 00 is 1, jump. 

If bit 00 is 0, load number 4000. 

Take logical product of LSH and 4000. 

If bit 11 of LSH is O, no action Bequineds 
If bit Jlof Loe ig 1, load. 3 ja. 

Make bit 11 of LSH equal to O. 

Store modified LSH. 

Unconditional jump. - 


If bit 00 of MSH is 1, load 4000. 

Is bit 11 of LSH 1 or 0? 

If 1, no action required, jump. 

If 0; load 4000. 

Make bit 11 of LSH 1 by adding 4000. 
Store modified LSH. | 


This last sequence of instructions has 
made the most significant bit (bit 11) 
of the LSH agree with the least sign- 
ificant bit (bit 00) of the MSH, thus 
removing the sign bit inserted in the 
two word product by the arithmetic unit. 


Load MSH. 

Left shift MSH 11 bits. This moves bit 
OO to the llth bit position for later 
removal (since it is no longer needed). 
Store modified MSH. 

Load the number 0777. 


Mask off upper three bits of LSH. 
Store lower nine bits, 
Load the number O400. 


Subtract lower nine bits of LSH. 
If result positive, no round off required. 
If negative, load number 7000. 


Mask off lower nine bits of LSH. 
Store modified LSH. 
Reload number 7000, 





SUBROUTINE KMOD (con't) 


Cell Mach Symbolic 
Code 
Re 3476 abd 
112 6102 n2f 
eo 5477 aod 
114 2200 eo ldf 
115 1000 - 
116 4076 add 
ly, 4076 std 
120 2200 f ldf 
ent 0777 
122 1077 lpd 
123 0110 sha 
124 4077 std 
125 2076 ldd 
126 0110 sha 
127 4076 std 
140 0407 ldn 
131 1076 lpd 
132 3077 add 
133 4167 sti 
134 2166 ldi 
135 6205 pjf 
136 2200 laf 
147 1462 
140 4167 adi 
141 4167 sti 
142 2073 i ldd 
143 6003 zjf 
144 2567 lei 
145 4167 sti 


pronde 
kadd 
at 


1462 
stoadd 
stoadd 
sign 
ind 
stoadd 
stoadd 
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Renarks 


Subtract LSH, to see if round-off correction 


will propagate to MSH. If not, jump. 
If so, add 1 to MSH. 
Load number 1000. 


Add 1000 to LSH for round-off correction. 
Store modified LSH. 


Since the least significant 9 bits of the 
LSH are discarded in the condensation of 
the two word product to a 12-bit word, a 
test is made to determine if these bits will | 
force a round-off correction to be made. 
If a round-off correction is to be made, 

a further test is required to determine 

if this correction will propagate to the 
MSH. Corrections are then made to the LSH 
and MSH if necessary. - 


NAPE FRAN Ona Ue 


eer eee. Tee eT Serr mee Le 


Load the constant 0777. 


Mask off upper 43 bits of MSH. 
Shift MSH left three bits. 
Store modified MSH. 

Load LSH. 

Shift LSH left three bits. 
Store. . 
Load the number 0007. 
Mask off upper nine bits of LSH. 
Add MSH to LSH to form K. 
Store K. 
Reload K' (unconverted). 

If positive, jump. 

If negative, load the number 1462. 


Add "2" to converted K. 

Store. 

Load sign of K', 

If positive, no action required. 
If negative, complement K. 
Store K. | 


The extraneous bits in the LSH and the 
MSH are set to zero and the remaining bits 
are shifted so that the LSH and the MSH 
can be added to form the 12-bit converted 
K. If 2000 was subtracted earlier, it is 
added back to the K in the form 1462. | 
The sign of K is made to agree with that 

of K', 





Cell 


146 
147 
150 
151 
152 
Noe 
154 
155 
156 
157 
160 
161 


Mach 
Code 
2070 
0701 
4070 
6010 
0402 
5067 
0402 
5064 
0402 
5066 
1071] 
7072 
0064 
0065 
0066 


' 0067 


0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0000 


. ind 


sinadd 
temp 
kadd 
stoadd 
index 
repadd 
link 
sign 
one 
outbuf 
leh 
msh 


SUBROUTINE KMOD (con't) 


Symbolic 
ldd index 
sbn Ol 
std index 
wD +5 
idm 02° 
rad  stoadd 
ldn 02 
rad sinadd 
ldn 02 
rad  kadd 
jpi repadd 
jpi - link 
equ 64 
equ > 65 
equ 66 
equ 67 
equ 70 
equ /l 
equ. 72 
equ 73 
equ. 74 
equ 15 
equ. 76 
equ 77 
end 
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Load number of (KTY)ts to be converted. 
Subtract l. 

Store. 

If zero, return to main progran. 

If not zero, load 0002. 

Advance storage address by 2. 


Advance sign address. 


Advance address of K' to be converted. 
Repeat subroutine for next K'. — 
Return to main program. 

Address of the sign of K!. 

Temporary storage cell. 

Address of K'., 

Storage address of K. 

Number of (K')'s to be converted. 
Address for subroutine repeat. - 

Link address to main program. 

Sign of K' presently under conversion. 
0651. 

Output location for K'. 

LSH of two word product. 

MSH of two word product. 





APPENDIX IV 


PROGRAM DIGIT 
IV-1. GENERAL DESCRIPTION. 

Program DIGIT is a CDC 160 digital computer program written for the 
control of the simulated hybrid oe. discussed in Chapter IV. A block 
diagram of the simulated system is shown in Figure 4-6. Specifically, 
DIGIT computes 

U(k) = +a, [X,(k)-- x, - fe(k)] + a,[X,(k) - 4] (IV-1) 
where X, (k); X, (k) » ft(k) afd r are digital representations (in the analog 
digital psuedo-number system) of sampled analog quantities received from 
the ADC and aj» a5 and r are manually-entered quantities. The quantity 
U(k) is transmitted to the DAC after computation is completed. All of the 
quantities in Equation IV-1 are further defined in Chapter IV. 

DIGIT uses subroutine KMOD to convert the quantities aj» a5 and ty 
from the normal octal representation of a decimal number to a form consis- 
tent with the analog-digital psuedo-number system generated by the A/D-D/A 
conversion process. AILI1 arithmeti¢ operations are then carried out in this 
number system using the CDC 168 arithmetic unit for multiplication and 
division. . 

For brevity, the following terminology is used in this discussion: 

(1) SUMI 2 CUS oa rt (k) 


(2) SUM2 = X, (k) -r. 


(3) P, = PRODI = +a, (SUM1). 


1 
(4) P, = PROD2 = +a,(SUM2). 


(5) SUM3 = P, aE P, = U(k). 


The following features are included in DIGIT: 
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(1) Overflow tests are performed on He P, and SUM3 and if over- 
flow occurs, automatic overflow correction or problem termina- 
tion is effected at the option of the user. Problem termina- 
tion includes the "resetting" of the analog computer and a 
visual indication of tthe quantity which overflowed. 

(2) A variable length time-delay loop is included permitting 
internal control of the sampling rate. This loop may be by- 
passed entirely if external timing control is desired. 

(3) The gperabeereset modes of the PACE TR-20 analog computer may 
be controlled settely by DIGIT. 

The various control options and data for subroutine KMOD are entered 
manually in low core cells. Execution of the program does not destroy this 
manually entered data, hence the program may be re-run without re-entry of 
data. 

DIGIT is written in a general format so that it may be used with any 
analog computer simulation of a second-order system where the control law 
consists of a linear combination of the states of the system. Any of the 
quantities in Equation IV-1 may be set to zero either by setting the parti- 
cular A/D channel at minus five volts for the sampled quantities or enter- 
ing a zero for the manually entered quantities. 

IV-2 USE OF PROGRAM DIGIT. 

A. LOAD program DIGIT and subroutine KMOD where desired in: the CDC 

160 computer memory. DIGIT requires 327 (octal) storage cells; KMOD re- 


quires 161 (octal) cells. In addition, low core cells 0000 thru 0013 and 


0040 thru 0077 are used by DIGIT. 
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B. ENTER the following data MANUALLY into the indicated locations: 


CELL NUMBER CONTENTS 
0000 Load address of KMOD. 
0001 0003 (number of conversions for KMOD). 
0002 0040 (memory address of r. after.conyer- 
sion iby) KMOD). 
0003 sign of r_ =}0000 if r_ is positive. 
O F ; : 
0001 if r, is negative. 
0004 |r | in x.xxx format (jr | SS). 
O O 
0005 \ Sign of a, = 0000 if a, is positive. 
0001 if a, is negative. 
0006 c Ja, in x.xxx format ( |a, | <<" 5)8 
0007 : Sign of a, =(0000 if a, is positive. 
2 : ; : 
0001 if a, is negative. 
0010 |a5| in x.xxx format (ja,| <1855)): 
0011 0000 if time delay loop is to be bypassed. 
0001 if time loop is to be used. 
0012 Coarse delay factor. See discussion below. 
0013 Fine delay factor. See discussion below. 


C. External Equipment. DIGIT requires CDC 160 access to the following 


external equipment: 


(1) CDC 168 arithmetic unit. 


(2) A/D-D/A conversion equipment with 


(a) 
(b) 
(c) 
(d) 
(e) 
(f) 


Xs available on channel #1, A/D. 


Xo available on channel #2, A/D. 

rt(k) available on channel #3, A/D. 

r available on channel #4, A/D. 

channel #1, D/A available for CDC 160 output. 

all input/output quantities appropriately biased so as 


to lie within zero to minus ten volts. A common bias is 


recommended. 65 
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(3) Remote control connections of the operate-reset switch of 
the PACE TR-20 analog computer (optional). 
D. RUN DIGIT from the initial load address. 
IV-3 OVERFLOW PROVISIONS. 
Product overflow detection is done by comparison of the product magni- 


tude a, x SUM and the variable magnitude SUM, | » provided that a 


i i 


is greater than unity. If the product magnitude is less than the variable 
magnitude, overflow has seeteeea. No test is made if a. is less than unity 
since overflow is not then possible. If overflow is detected, the computer 
will either error halt (providing a visual indication as to which product 
overflowed) or will substitute the maximum possible value (4000 or 3777) 
for the intended product depending upon the signs of a; and SUM, . If the 
error halt option is used, tile analog computer is reset to the initial 
conditions. 

Summation overflow detection is accomplished by detection of the sign 
change of SUM3 as overflow occurs, provided that both arguments of the sum 
agree in sign. If the arguments are of opposite sign, no overflow is possi- 
ble. The same provisions concerning error halt or overflow correction are 
provided as for product overflow. 

If the error halt option is selected, an indication of the quantity 
causing the overflow is displayed in the A-register of the 160 computer. 
The error halt may be bypassed by manual entry of the appropriate error by- 
pass code into the A-register (after clearing of the registers) and re-run- 
ning of the program. This will provide automatic owerflow correction. The 


error halt indications and error bypass codes are tabulated in Table IV-1l. 
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Table IV-1. Overflow Error Halt and Bypass Table. 


IV-4. INTERNAL TIMING CONTROL. 

The variable-length time delay loop in program DIGIT is:shown in 
flow chart form in Figure IV-1. Delay is effected by loading the fine delay 
factor into the arithmetic register of the CDC 160 computer and counting 
this number down to zero by successive subtraction of the number one. This 
count-down procedure is repeated as many times as is indicated by the coarse 
delay factor. A total of 52.4 milliseconds is required to count down the 
number 7777 to zero. A suggested method for accomplishing a time delay of 
T milliseconds is: 

(a) Set the fine delay factor to 7777. 

(b) Set the coarse delay factor to the first octal digit exceed- 
ing T/52.4. (T in milliseconds). This should yield a time 
delay slightly over T milliseconds in length. Using an ex- 
ternal time-measuring device (oscilloscope or Brush recorder), 
adjust the fine delay factor downward to yield the desired time 


delay. 
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Example: To achieve a one second sample interval for program DIGIT, the 


coarse delay factor was set to 23, and the fine delay factor was adjusted 


8 
to 7600... 


8 
IV-5 EXTERNAL TIMING CONTROL. 

Insertion of zero in cell 0011 will bypass the internal timing loop of 
program DIGIT. External ctifhing “CREE may be accomplished by the connec- 
tion of an external clock to the "input disable" jack of the ADC. This jack 
is connected via an AND-gate to the “input ready'' line of the 160 computer 
input cable, The 160 computer will delay execution of the input instruc- 
tion until the “input ready" line drops to -13 volts. Thus, if the “input 
disable" is held at ground level by the external clock but periodically 
dropped to -13 volts Bt chewcoan iets frequency, sample delay is accomplished. 
Since input control is accomplished on the "input ready" line vice the “input 
request'' line, the first sample that is inputted to the 160 must be dis- 
carded since this sample is stale. Also, the “input disable" must be held 


at -13 volts long enough for four channels to be sampled each period (about 


450 microseconds). 
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csdel = coarse delay 
factor loaded in cell 
0012. 


load fndel 
fndel = fine delay 
fndel = fndel - 1 factor loaded in 


om cell 0013 
fndel =0 > 
Ss 4 
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del = del -l 


Yes 


Return to 
main program 


Figure IV-l. Program DIGIT Variable-Length Time Delay Loop 
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IV-6. DIGIT Flow Chart and Program Listing. 

Figure IV-2 is a flow chart presentation of Program DIGIT. Pages 
72 ‘thru 77° contain a program listing of DIGIT as compiled ee the 
OSAP compiler for the CDC 160 digital computer. Both machine language and 
mnemonic listings are presented. 

It should be noted that all arithmetic operations within the sample 
loop are carried out using the analog-digital psuedo-number system as 
explained in Chapter 3. This method requires a division by (0631), after 


each product is formed. 


4 
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Figure IV-2. DIGIT Flow Chart 
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Cell Mach Symbolic Remarks 
Code ; 

0100 4047 std errflg Store error stop code. 

0101 -0101 sha Ol Oall subroutine KMOD. 

0102 7000 jpi sradd | 

0103 2200 ldf 

0104 0631 631 

0105 4062 std  onea Store 0651, 

0106 4065 std oneb 

0107 O10] sha Ol Generate and store 

0110 0603 adn 03 loop address. 

0111 4075 std loop 

0112 . 7500 exf 

0113 1401 1401 

0114 7265 inp chil Input Xj. 

0115 0077 Ah a 

0116 7500 exf 

0117 1402 1402 

0120 7262 inp ch2 Input X5. 

0121 0100 ~ 100 

0122 7500 exf 

0123 1403 1403 Input ramp = 

0124 7257 inp ch3 rdot x t(k). 

0125 0075 iS 

0126 7500 exf 

0127 1404 1404 : 

0130, 7254 inp  ch4 Input rdot. 

014] 0074 74 

0132 2076 ldd@ xl 

0133 © 3474 sbd ramp Compute suml e 

0134 3440 sbd step X)} = ramp - step. 

0135 4043 std suml 

0136 2077 ldd = x2 Compute sum2 = 

0137 3473 abd  rdot Xo - rdot. 

0140 ©4045 std sum2 

0141 7500 exf 

0142 3301 3301 

0143 = =7330 out al 

0144 0044 4h 

0145 7230 inp pl Compute 

0146 0062 62 P) = a) x sunl, 

0147 7500 ext 

0150 3300 3300 ay = -1 for 

0151 7324 out pl optimum control. 

0152 0063 63 

0153 7224 inp ql 

0154 0070 70 


PROGRAM DIGIT 
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PROGRAM DIGIT (con't) 








Cell Mach Symbolic Remarks 
Code 
O9S> . (500 exf 
Gi5ome 5301 3301 
0157 135 out a2 
0160 0046 46 
0161 7215 inp p2 Compute 
0162 0065 65 Po = ap xX sum. 
0163 7500 exf \ | 
0164 3300 3300 ao = -1.5 for 
OVER 7321 out p2 optimum control. 
0166 0066 66 
0167 7211 inp q2 
0170 0072 12 
0171 6114 nzf- 64 
0172 6013 2jf a 
ese eeeeeee ee se® & = = Sr! a ae ne ee ee ee ee ee Se 
0173 oo42 al 42 
0174 ° 0044 a2 hd 
0175 0060 pi 60 
0176 0063 p2 63 
0177 0066 ql ~ 66 
0200 0070 q2 . 10 Input-output table. 
0201 0076 chi 76 
0202 OO77 ch2 TL 
0203 0074 ch3 74 
0204 0073 ch4 i, 
0205 2043 a ldd = suml . . 
0206 6205 Dita OD Determine sign and 
0207 2443 led = uml magnitude of suml. 
0210 4057 std magsl 
0211 odo] ldn 01 
0212 6103 nzf it 
0213 4057 »b std magsl 
0214 0400 ldn 00 
0215 4056 +t std sinsl 
0216 2045 ldd sum2 
0217 6205 Dit ec 
0220 2445 led sum2 Determine sign and 
0221 4055 std mags2 magnitude of sum2. 
0222 O40] lida Ol 
0223 6103 nzf ou 
0224 4055 ¢ std mags2 
0225 0400 lIdn 00 
0226 4054 u std sins2 
0227 2066 ldd  prodl 
0230 6205 pjf | 
0231 2466 led = prodl Determine sign and 
0232 ° 4053 std magpl magnitude of P). 
0233 Oo401 ldn (0) 
0234 6103 nzfoooyv 
0235 4053 4 std magpl 
0236 0400 ldn 00 
0237 4052 v std  sinpl 


5) ae 








Pelaisi. 
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PROGRAM DIGIT (con't) 


Cell Mach Symbolic Remarks 
, Code 
0240 2070 Idd prod2 
0241 6205 pif ie 
0242 2470 led prod2 Determine gign and 
0245 4051 etd magp2 magnitude of Po. 
(0244 401 ldn Ol 
0245 6103 n2f ow 
0246 4051 e std magp2 
0247 o400 ldn 00 
0250 4050 w std sinp2 
0251 2006 ldd alu 
0252 -1252 lpf const 
0253 6033. Bite fT 
0254 205% ldd  magpl P) overflow test. 
0255 3457 sbd magsl 
0256 6230 pje f 
0257 6027 ne | 
0260 2047 ldd errflg 
0261 0201 ipn of 
0262 6105 nzf sg 
0263 7500 exf 
0264 2401 2401 Error stop 1 if 
0265 0401 ldn 01 P) overflows. 
0266 0000 err 
0267 2056 g ldd  sinsl 
0270 6004 2if oh 
0271 2005 ldd  sinal 
0272 6004 zjf ce 
0275 G07 nzf dd 
0274 2005 h ldd = sina] 
0275 6005 2jf dd 
0276 O401 ~ ce ldme Oi P| overflow correction. 
0277 4052 std  sinpl 
0200 2226 ldf nmax 
0301 6104 n2f aa 
0302 o400 dd ldn 00 
0303 4052 std  sinpl 
0304 2221 ldf pmax 
0305 4066 aa std prodl 
0306 2010 f ldd a2u 
O507¢—5 12> lpf const 
0310 6036 Za 
OF1l8” Zeal ldd magp2 Po overflow test. 
0312 3455 sbd mags2 
0315 6233 pe 5 
0314 6032 zafe 5 
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PROGRAM DIGIT (con't) 


é 


Cell Mach Symbolic Remarks 





Code 
0315 2047 ldd orrfig . 
0316 0202 lpn 02 
0517 6110 nef ik 
0720 F000 exf 
0321 2401 ‘2401 Error stop 2 if 
0322 0402 ldn 02 Po overflows. 
0323 0000 | err 
0324 7000 const 7000 
0525 3/// pmax oii Table of constants. 
0326 4000 nmax 4000 
03527 2054 &k ldd  sins2 
0330 6004 mit 1 
0331 2007 ldd  sina2 
0332 6004 Zc ee 
0355 «6107 -~ nef ee - 
0334 2007 #1 ldd  sina2 Po overflow correction. 
03355 6005 zjf ee 
0336 o401 ff ldn 01 
0337 4050 std sinp2 
0340 ©2312 ldb  nmax 
0341 6104 nzf bb 
0342 0400 ee ldn 00 
0343 4050 std sinp2 
0344 2317 ldb pmax 
0345 4070 bb std prod2 
03246 2066 j ldd = prodl 
0347 3070 add prod2 Compute 
0350 4072 std sum sum) = Py # Po. 
0351 2052 ldd = sinpl 
0352 43450 sbd sinp2 
0353 6125 nzf mm 
0354 2072 ldd sum Determine sign of 
0355 “6205 pat iin sum} and perform 
0456 o401 ldn 01 overflow test. 
0357 6102 nzf oo 
0360 O400 n lan 00 
0361 3452 o sbd sinpl 
0362 6016 zjf om 
0363 2047 ldd errflg 
0364 0204 lpn 04 Error stop 34 if 
0365: 6105 nzf sp sum? overflows. 
0366 7500 exf 
0367 2401 2401 
0370 0403 ldn 03 
0371 0000 err 
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PROGRAM DIGIT (con't) 


Cell Mach Symbolic 
Code 
0372 2052 _ »p ldd 
0373 6003 Za, 
0374 2346 ldb 
O5/75  OfGe2 n2r 
0376 2351 q ldb 
0377. 4072 Fr std 
0400 7500 am exf 
0401 2411 
0402 7325 out 
0403 0073 
0404 2011 ldd 
0405 6013 2jf 
0406 2012 ldd 
0407 4046 std 
0410 20143 ol _ ddd 
0411 0701 vit sbn 
0412 6501 nzb 
0413 2046 ldd 
0414 0701 sbn 
0415 4046 std 
0416 6506 nzb 
0417 7075 8 jpi 
0420 7500 zz exf 
0421 1401 
0422 7204 inp 
0423 odode 
0424 6505 nzb 
0425 6406 zjb 
0426 0041 dump 
0427 0072 sum 
0000 sradd equ 
0005 sinal equ 
0006 alu equ 
0007 sina2 equ 
0010 a2u equ 
OOll delflg equ 
0012 ~ =ecsdel equ 
0013 fndel equ 
0040 step equ 
0043 = suml equ 
0045 sum2 equ 
0046 =del equ 
0047 errflg equ 
0050 sinp2 equ 
0051 #£magpe equ 
0052 ~=s siinpl equ 
0055 #£magpl equ 


Remarks 
sinpl 
q 
nmax 
r 


Sum3Z overflow correction. 


pmax 


sum? 


Delay loop. 


Dump lst sample if external 
timing used. 


Symbol-cell equivalence 
table. 
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PROGRAM DIGIT (con't) 





Cell Mach Symbolic Remarks 
Code 
0054 sins2 equ 5] | 
0055 #£mags2 equ 55 
0056 ~—s sins equ 56 
0057 #£magsl equ 57 
0062 onea equ. 62 
0065 oneb equ 65 
0066 —s prod equ 66 
0070 ~=prod2 equ. 70 
0072 sum) equ {2 
0073 =rdot equ 7/3 
0074 ramp equ 74 
0075 loop equ 75 
0076 Xl equ 76 
0077 X2 equ 77 
0000 end 
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